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Abstract 

Recent works on MIMO interference channels have shown that interference alignment can 
significantly increase the achievable degrees of freedom (DoF) of the network. However, most of 
these works have assumed a fully connected interference graph. In this paper, we investigate how the 
O ' partial connectivity can be exploited to enhance system performance in MIMO interference networks. 

We propose a novel interference mitigation scheme which introduces constraints for the signal 
J> . subspaces of the precoders and decorrelators to mitigate "many" interference nulling constraints 

QQ ' at a cost of "little" freedoms in precoder and decorrelator design so as to extend the feasibility 

' region of the interference alignment scheme. Our analysis shows that the proposed algorithm can 

l/^ I significantly increase system DoF in symmetric partially connected MIMO interference networks. 

o 

We also compare the performance of the proposed scheme with various baselines and show via 



simulations that the proposed algorithms could achieve significant gain in the system performance 
of randomly connected interference networks. 

I. Introduction 

Recently, there is an intense research interest in the area of interference channels and the 
interference mitigation techniques. Interference alignment (lA) was proposed in [21, ^ to 
reduce the effect of multi-user interference and is extended to deal with interference in MIMO 
X-channels O and K pairs interference channels [0). The key idea of lA is to reduce the 
dimension of the aggregated interference by aligning interference from different transmitters 
into a lower dimension subspace at each receiver. Using infinite dimension extension on 
the time dimension (time selective fading), it is shown that the lA could achieve the optimal 
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Degrees-of-Freedom (DoF) of in fT-pair MIMO ergodic interference channels [HJ with N 
antennas at each node. In |[3, the authors proposed the concept of ergodic alignment, which 
also utilizes symbol extension exploiting time- selective fading of interference channels. 

One important challenge of lA scheme is the feasibility condition. For instance, the lA 
schemes in BU requires 0{{KNY^^^^) dimensions of signal space to achieve the total 
DoF. To avoid such huge dimensions of signal space, some researchers have studied lA 
designs for quasi-static MIMO interference channels. With limited signal space dimensions, 
the achievable DoF of each transmitter-receiver pair in MIMO interference channels is upper 
bounded by (where K is the number of transmitter-receiver pairs, Nt, are the 

number of antennas at each transmitter and receiver, respectively) BH. Unlike the time- 
selective or frequency-selective MIMO interference channels, total DoF of quasi-static MIMO 
interference channel does not scale with K. Furthermore, it is quite challenging to design 
precoders and decorrelators that satisfy the lA requirements in limited dimension MIMO 
interference channels due to the feasibility problems [|6]. In an iterative precoders and 
decorrelators design based on alternating optimization is proposed for quasi-static MIMO 
interference channels. In flSl, flU, some constructive methods to design precoders and decor- 
relators are proposed, but these schemes can only achieve 1 DoF per transmitter. 

In fact, the technical challenge on the feasibility issue in limited dimension MIMO interfer- 
ence channels is highly related to the full connectivity in the interference graph. In practice, 
the interference channels are usually partially connected due to path loss, shadowing as well 
as spatial correlation. Most of the existing literatures have assumed fully connected MIMO 
interference channels such as equal path loss and spatially uncorrelated MIMO channels. 
Intuitively, partial connectivity may contribute to limiting the aggregate interference and this 
may translate into DoF gains in the system. In this paper, we are interested to study the 
potential benefit of partially connectivity in a i^'-pair MIMO interference network with quasi 
static fading. There are several important technical challenges involved. 

• How to exploit partial connections in interference mitigation? Traditionally, it is 
well-known that partial connection (due to path loss, shadowing or spatial correlation) 
is detrimental to point-to-point MIMO performance [[TOl . [fTTTl because it reduces the 
number of spatial channels in the MIMO link. However, in MIMO interference chan- 
nels, partial connection may also reduce the dimension of the undesired signals (the 
interference), leading to possible performance improvement. In other words, we can 
potentially design precoders to exploit the partial connection property and reduce the 
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interference dimensions to other users. 
• Achievable DoF for partially connected MIMO interference channels. In order to 
obtain insights on the potential benefits or degradations of partial connections in MIMO 
interference channels, one would be interested in deriving DoF bounds. Existing DoF 
analysis of interference channels 0), [|6l assumed i.i.d. MIMO fading channels (fully 
connected MIMO interference channels) and it is interesting to find out how the partial 
connection parameters such as node density and spatial correlation level affect total DoF 
of the system. 

In this paper, we propose a novel two-stage dynamic interference mitigation scheme to 
exploit the potential benefit of general partial connections in limited dimension MIMO inter- 
ference channels so as to improve the network total DoF. The proposed dynamic interference 
mitigation solution has two stages. The first stage determines the stream assignment and the 
subspace constraints for the precoders and the decorrelators based on the partially connected 
topology such as the path loss, shadowing and spatial correlation. The second stage determines 
the precoders and the decorrelators (based on the instantaneous channel state information) 
over the subspaces obtained from the first stage. Based on the proposed dynamic interference 
mitigation scheme, we shall derive an achievable DoF bound of a symmetric interference 
network and show that the DoF in partially connected MIMO interference channels can exceed 
the well-known DoF results of for i.i.d. MIMO interference channels. Furthermore, we 

shall discuss how the DoF gain is affected by the partial connectivity (path loss and spatial 
correlation) in the system. Finally, we shall compare the performance of the proposed scheme 
with various conventional baselines, and it can be observed that the proposed scheme offers 
significant performance gain over a wide range of system operating regimes. 



A. General K-pair Partially Connected Quasi-static MIMO Interference Channels 

We consider a MIMO system with K transmitter (Tx) and receiver (Rx) pairs. Each 
transmitter and each receiver has Nt and A^^ antennas, respectively. Denote the channel fading 
coefficients from the Tx n to the Rx m as H^n ^ C^''^^*. Let dm {< min(A't, A^,.)) be the 
number of data streams (DoF) transmitted by Tx-Rx pair m. The received signal G C^™ 
at Rx m is given by: 



II. System Model 




(1) 
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where G C^™ is the encoded information symbol for Rx m, G C'^™^^'^ is the 
decorrelator of Rx m, and G C^'^"'" is the transmit precoding matrix at the Tx m. 
z G C^"^^^ is the white Gaussian noise with unit variance. The transmit power at the Tx n 
is E(||V„x„|p) = P„. The channel connectivity of the fC-pair interference channels {Hm„} 
is specified by the following model. 

Assumption 2.1 (General Partially Connected Model): The elements of the channel states 
matrices {Hmn £ C^"^^^'}, n,m E {1,2, K} are random variables following certain 
distribution and have the following properties: 

• Independence: Random matrices n,m E {1,2, K} are mutually independent. 

• Partial Connectivity at the Transmitter Side: Define the Tx side partial connectivity 
as the null space of Hmn, e-g: A/'(Hm„) = {v G C^'^^ : H„j„v = 0}. 

• Partial Connectivity at the Receiver Side: Define the Rx side partial connectivity as 
the "transposed" null space of U^n, e.g: A/'^(H^„) = {u G C^^^' : H^„u^ = 0} = 



As a result, {J\f{ii„in)} and {A/'^(H,mn)} m,n E {1,2, K} are the connection topo- 
logical parameters of the general partially connected model. This is a general model as no 
specific structure is imposed on {MiYlmn)} and {A/'^(Hmn)}- To help readers get some 
concrete understanding on the physical scenarios when we have partial connectivity, we shall 
illustrate some examples of partially connected interference channels in the next section. Note 
that the partial connectivity model imposed in Assumption 12.11 contains, but is not limited 
to, these examples. 

B. Example Scenarios of Partially Connected Systems 

1) Path Loss and Shadowing: In practice, different Txs may contribute differently to 
the aggregate interference due to the heterogeneous path loss and shadowing effects. For 
example, in a /C-pair MIMO interference network with Nt = Nj., suppose some Txs and Rxs 
are far away from each other when the difference between their node indices \n — m\ > L, 
the path loss and shadowing from Tx n to Rx m, is 60 dB higher than that of the direct 
link (from Tx n to Rx tt,, n G {1,2, K}) and hence, effectively, we have H^n = 0, 
V|n — m| > L. This corresponds to a partially connected MIMO interference channel (induced 
by path loss and shadowing effects) with the connection topology given by A/'(Hmn) = 



{u G C 



: uH, 



mn 



0}. 
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2) Unequal Transmit and Receive Antennas: When Nt ^ Nr, there is a rank | A't — A^^l 
null space on the side with more antennas. For example, when A^^ = 4, A'^ = 2, denote 

Hmn = ^ where h.mn{p) are 1 X 4 vectors. Hence, this corresponds to a partially 

hmn (2) 

connected MIMO mterference channel (induced by non-square fading matrices) with the 
connection topology given by: A/'(H„„) = (span(h^„(l), h^„(2)))-L, A/'^(H^„) = {0}, 
Vm,n G {1,2,...,K}. 

3) Spatial Correlation: As shown in [fT2ll . lfT3]| . in practice, local scattering effect causes 
significant spatial correlation in MISO channels. To explore the similar effect in MIMO 
channels and get a direct association between the MIMO fading channel correlation and the 
physical scattering environment, we shall introduce the virtual angular domain representation 
for MIMO channels [fTT]|. Specifically, the MIMO fading channels in the antenna domain 

and in the angular domain H'^-} have a one-one correspondence given by: 

H^^ = A^Hr A^ (2) 

where: At = [ejv,(0), e^v, (^-^^ ...,ejv, (3) 
A„ = K,(0),e„,(-l).^..e.,(^)l (4) 
e„(uj) = 1 [i,e-i2-(")_e-''^-<2")...e-J2-«"-")")]'- (5) 



Given a local scattering environment with the parameter effective scattering radius S as 
illustrated in Fig. [L the MIMO fading matrix = {h'^i{p,q)} p G {1,2, ...,iV,.}, q e 

{1,2,..., Nt} in the angular domain has the following property: h^fl{p, q) = with probability 
1 if and only if (Please refer to Appendix [A]) 



sin 6 q 



mod 1 > — , e [9mn-Fa + Fa{S, d^n)] (6) 



where Fa{S, d„ 



arcsin when: S < dmn 
TT when: S > dmn 



dmn £ (~^) is the direction from Tx n to Rx m, and d^n is the distance between the 
two nodes. This spatially correlated MIMO model is a special case of the general partially 
connected MIMO interference model. For example, in Fig. \T\ suppose A^^ = A^^, = 4 
and denote H^^ = [h^„(l), h„„(2), h„„(3), h„„(4)]. Suppose due to spatial correlation, 
hmn(l) = hm„(4) = and hmn(2), hmn(3) are randomly generated C^^^ vectors. AfCHmn) = 
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span 



1 

i 




U 























1 



V 

Vn,m G {1,2,...,K}. 



and A^^CH, 



(span {{ARhmn{2)f, (A^h„„(3))^))- 



III. Algorithm Description 

In this section, we shall propose a novel dynamic interference mitigation scheme to exploit 
the topological advantage due to partial connectivity. The proposed scheme is also backward 
compatible with existing lA designs when the topology is fully connected. The algorithm 
dynamically determines the data stream assignment © = {di, d2, d^}, dn G {0, 1, rf™'*''} 
and the associated precoders V„ G C^*^'^" and decorrelators U„ G C''" x A^,,, where d™^^ and 
dn are the number of the data streams claimed by and assigned to Tx-Rx pair n, respectively, 
n G {1, 2, K}, such that: 



rank(U„H„„V, 



dr. 



0,Vn^mG {1,2, ...,K} 



(7) 
(8) 



Most of the existing works on lA have assumed a fully connected interference topology, e.g. 
{Hm„} are all full rank. However, as we have illustrated in Section ITl-Bl MIMO interference 
channels are usually partially connected due to various physical reasons. As far as we are 
aware, no existing schemes can be extended easily to exploit the potential benefit of partially 
connectivity in MIMO interference networks. 



A. Motivations: A Dynamic Interference Mitigation Scheme for General Partially Connected 
MIMO Interference Channel 

1 ) The potential benefit of partial connectivity: We shall first illustrate the potential 
benefit of partial connectivity by a simple example. 

Consider a 2 x 2, 5-pair interference network. Each Tx-Rx pair attempts to transmit 1 data 
stream (i.e. d™^"" = 1, G {1,2,3,4,5}). If the network is fully connected (all channel 
matrices are rank 2), the freedoms in each precoder G C^^^ and decorrelator G C^^^ 
is given by dim(^(l,2)) = 1, where the Grassmannian [[T6l . [[T71 Q{x,y) denotes the set 
of all x-dimensional subspaces in C^, and the number of constraints induced by each cross 
link (UmHmnVn = 0, H ^ m) IS 1. If we assign data streams to k Tx-Rx pairs, there are in 
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total 2k freedoms in the precoders and the decorrelators and k{k — 1) interference alignment 
constraints. Hence, from the lA feasibility condition [6|, we have k{k — 1) < 2k ^ k < 3. 
In other words, the achievable DoF is upper bounded by 3. 

Now suppose the network is partially connected such that the channel matrices of all cross 
links are rank 1 with null spaces given by the red arrows in Fig.[3l Then we have lA constraints 
(|7]), dH]) are satisfied under the following policy: Assign data streams to Tx-Rx pairs 



({1,2, ...,K}={1,2,4,5}), with precoders vi = ^ 
, and decorrelators Ui 

H42V2)^, Us = 



1 




1 












-1 










V5 



U4 





1 

-1 

1 



-1 

1 



-1 

1 



V2 

2 



,2,4,5 
1 

1 



-1 

1 



H21V1 



H54V4)^. Since 4 Tx-Rx pairs can have data 



streams simultaneously, 1 extra DoF is achieved compared to the fully connected case. 

The example above illustrates how partial connectivity can contribute to network perfor- 
mance gains. However, as we shall explain below, there are various technical challenges to 
exploit the benefit of partial connectivity for general scenarios. 

2) The difficulty in exploiting the benefit of partial connection: 
• Interference overlapping versus interference nulling: Classical interference alignment 
schemes |[4l reduce interference dimension by "overlapping" the interferences from dif- 
ferent Txs. However, when partial connectivity is considered, "overlapping" interferences 
is no longer the only method to reduce interference dimension. Part of the interference 

can also be eliminated by utilizing the null spaces of the channel states. For instance, 

1 

, we eliminate the interference 



V2 
2 



-1 



in the previous example, by setting vi 

from Tx 1 to Rx 4 and 5 by utilizing the null space of H41 and H51. In practice, we 
should dynamically combine these two approaches in order to better exploit the nulling 
opportunities as well as alignment opportunities in a partially connected interference 
network. However, a combined design may depend heavily on the specific realization of 
the connection topology and a combination criteria that can work for general connection 
topologies is not yet clear. 
• Freedoms versus constraints: Another perspective of the technical challenges is on the 
feasibility conditions in quasi-static MIMO interference networks. In [O, a symmetrical 
MIMO interference system is feasible if and only if the number of freedoms in transceiver 
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design is no less than the number of independent constraints induced by interference 
alignment requirements. In order to reduce the number of the independent constraints in 
a partially connected network, we need to restrict the precoders and decorrelators to a 
lower rank subspace. On the other hand, restricting precoders and decorrelators to a lower 
rank subspace reduces the number of the freedoms in the precoders and decorrelators. 
The conflict between freedoms in transceiver design and independent constraints makes 
the subspace selection very challenging. 

• Exponential complexity in checking lA feasibility conditions: As revealed in [6|, 
checking the lA feasibility condition requires comparison of freedoms versus constraints 
for every possible combination of the interference nulling constraints. This process 
involves 2^(^-i) - 1 comparisons, where K is the number of Tx-Rx pairs. Such a 
complexity is intolerable in practice. Hence, a low complexity algorithm for checking 
the feasibility condition on a real-time basis is needed. 

B. Dynamic Interference Mitigation Scheme for a 5 -Pair 2x2 Partially Connected Interfer- 
ence Network 

The following observation is the key insight of the proposed algorithm: 
Observation: In a partially connected MIMO interference network, by properly restricting 
precoders and decorrelators U„ to a lower rank subspace, we can eliminate "many" 
independent constraints at a cost of only a few "free variables " and hence extend the lA 
feasibility region. 

In this section, we shall use the example of a 5-pair 2x2 partially connected interference 
network described in Section llTl-All to illustrate the main ideas of the proposed scheme. 

• Step 1 Initialization: Note that all the 5 direct links have sufficient rank (rank(H„„) = 
2 > 1, n G {1, 2..., 5}), initialize D = {1, 1, 1, 1, 1} to see if the network is feasible with 
all Tx-Rx pairs active. 

After the initialization, the number of the freedoms in transceiver design and the number 
of the lA constraints is illustrated in Fig. UK- The numbers in red and blue denote the 
freedoms in the corresponding decorrelators and precoders , respectively, and the numbers 
in purple denote the number of the lA constraints. 

• Step 2 Find out the common subspaces in partial connectivity: As indicated by 
the"i/" signs in Fig. [3l from Tx 1 to Rx 4 and 5, there is a one dimensional common 
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subspace in partial connectivity state: A/'(H4i) = A/'(H5i) = span . Similarly, 

-1 

from Tx 2 to Rx 1 and 5, from Tx 4 to Rx 1 and 2, from Tx 5 to Rx 2 and 4, there 
are common null spaces. 

Step 3 Select subspaces to reduce the number of lA constraints: Set vi = 



V2 
2 

_^ 
2 

we have H„iVi = 0, m G {4, 5}. Hence, as indicated by the highlight parts in Fig. |43, 



we reduce 2 constraints at a cost of 1 freedom. Similarly, as indicated by the highlight 
parts in Fig. |4]I!, since Tx 2,4,5 each has two cross links with overlapping null spaces, 
we can reduce 2 constraints at a cost of 1 freedom by setting the precoder vectors to 
be the basis vectors of the corresponding null spaces. 

Step 4 Check the feasibiUty conditions: The number of the remaining freedoms and 
constraints after step 3 is illustrated in Fig. |4}r. Randomly assign the constraints to the 
corresponding Txs or Rxs (as indicated by the color of the numbers, deep-blue and 
deep-red indicate the number is assigned to the Tx or the Rx, respectively), we get 
Fig. [5Kl. In this figure, some nodes are "overloaded" in the sense that the number of 
freedoms at this node minus the number of constraints assigned to this node is negative 
(highlighted using yellow) while some nodes still have extra freedoms (highlighted using 
green). Hence, as illustrated in Fig. |5}\2, we can reassign the constraints so that we can 
have less overloaded nodes (the changes are highlighted using red boxes). However, in 
Fig. [5K2, there are still some overloaded nodes while no nodes have extra freedoms, 
hence the network is not feasible. 

Step 5 Turn off the most "constraint demanding" stream: As illustrated in Fig. |5K2, 
if we remove Tx-Rx pair 1 from the active set {1, 2, K}, we reduce 1 freedom (the 
freedom in vi) and 4 constraints (link Tx 1 to Rx 2,3 and Tx 3,5 to Rx 1). Hence, the 
freedom-constraint gain by removing Tx-Rx pair 1 is (—1) — (—4) = 3. Similarly, the 
freedom-constraint gains by removing Tx-Rx pair 2,4,5 are also 3, whereas the gain by 
removing Tx-Rx pair 3 is (—2) — (—8) = 6. Since 6 > 3, we remove the stream of 
Tx-Rx pair 3, i.e. let D = {1, 1, 0, 1, 1} and return to Step 3. 

Repeat Step 3~4: After removing the stream of Tx-Rx pair 3, repeat Step 3~4 and the 
number of the remaining freedoms and constraints after subspace selection is illustrated 
in Fig. [5^1. As revealed by Fig. [5^2, the lA constraints can be assigned properly so 
that no node is overloaded, the network is feasible. Continue to Step 6. 
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• Step 6 Precoder and decorrelator determination: Based on the results given in Steps 
1~5, use the minimum interference leakage iteration []3 to determine the precoders and 
decorrelators and we get the policy proposed in Section IIII-All 
Remark 3.1 (Low Complexity lA Feasibility Checking in Step 4): To avoid the exponen- 
tial complexity in lA feasibility checking, we have proposed a low complexity method, namely 
the freedom-constraint assignment with worst case complexity 0{K^) only. Furthermore, we 
shall formally prove in Appendix |D] that this method is indeed a necessary and sufficient 
condition of the lA feasibility conditions (fTTl) in the general case. ■ 

C. Dynamic Interference Mitigation Scheme - General Case 

Inspired by the example above, we shall propose a two-stage algorithm, namely stream 
assignment and subspace determination stage and precoder / decorrelator determination 
stage. The first stage algorithm (which corresponds to Step 1~5 in the example illustrated 
above) determines the stream assignment pattern D = {di, c^a'} and the subspaces 
for the precoders S* and decorrelators for n, m G {1,2, K} based on the partial 
connectivity state {A/'(Hm„), A/'^(H„„)}, m,n E {1,2, ...,K}. Based on the outputs of the 
stage I algorithm and the channel state, the stage II algorithm (which corresponds to Step 
6 in the example illustrated above) determines the precoder s/decorrelators. In the next two 
subsections, we shall elaborate the details of the stage I and stage II algorithms, respectively. 

1) Stage I: Stream Assignment and Subspace Determination: Suppose the row vectors in 
precoder V„ and the column vectors in decorrelator Um are constrained to the linear spaces 
§^ C C^*^i and C Ci^^% respectively. Denote S"* = = |§[„|, where |X| denotes 

the rank of linear space X. Then (|7]) and ^ can be rewritten as: 

rank(U„H„„V„) = rank(U'„S;H„„S^V:,) = (9) 

UmHmnV„ = U^SJ^HmnS^V^ = (10) 

where the column vectors of the Nt x 5** matrix and the S^^ x A^^ rnatrix SJ'^ span the 
spaces and SJ^^, respectively. Note that the constrained precoder V'^ and the constrained 
decorrelator U'^ are S"* x dn and dm x matrices, respectively. 

Before we elaborate the details of the stage I processing, we shall first define the notion 
of a proper MIMO interference system. 
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Definition 3.1 (Proper MIMO Interference Systems for the General Case): The MEMO in- 
terference system is proper for interference alignment if: 



VGT,Gij C {1,2, G = (G/j X GT)\{(n,n) ,n G G/j n G^}, where "x" denotes 



Remark 3.2 (Physical Meaning of Definition liJI) ; As proved in Appendix |Bl the left hand 
side of (fTTl) represents the total number of interference alignment constraints of the links from 
a Tx in set G„ to a Rx in set G^, and the first and second term on the right hand side of (fTTT i 
represent the sum of the free variables in V'^, n G G„ and U'„, m E Gm, respectively. Hence, 
(fTTl) means that for any subset of Tx-Rx combination G„©Gm C {1,2,..., iC}©{l, 2, K}, 
the number of constraints is no more than the number of free variables. ■ 

The main steps of the Stage I processing algorithm for the general K-pair partially 
connected interference network is illustrated below. Steps 1~5 below corresponds to the 
Steps 1~5 for the 5-Pair example. 

Stage 1 (Stream Assignment and Subspace Determination): 

• Stepl Initialization: Initialize the number of stream assigned to each Tx-Rx pair to be 
the minimum of the rank of the direct link and the number of streams claimed by this 
Tx-Rx pair, i.e. dn = min(rank (H„„), Vn G {1, 2, K}. 

• Step2 Calculate the common null spaced] : For every Tx-n, n G {1, 2, K}, calculate 
the common subspaces of the null spaces of the cross links from this Tx within the 
effective subspace of the direct link , i.e. A/'({Hm„ : m G Ksub}) = (nmeK,„bA/'(Hm„))n 
(A/'(H„„))^, C {1,2, ...,K} as follows: 

- Denote K*„ = {m : m G {1, 2, K},^^{Umn) ^ C^*^i(i.e. ^ 0)}, Initialize 
Ar(0) = (Ar(H„„))^, Ar({H„„}) = Ar(H„„) n (A/'(H„„))^, Vm g K, and subset 
cardinality parameter C = 2. 

- For every Kgub ^ IK^ with \Ksub\ = C, if all the subsets of IKs^ib with cardinality (C— 
1) are not {0}, calculate Af{{B.mn ■ m G K^^;,}) = A/'({H„„ : m G ]Ks„fe\{m'}} n 

'The worst case complexity of this step is 0{K2^ ~^), where K' — max„(max(|K^|, |KJ^|)). In practice, due to path 
loss, K' usually does not scale with K. Hence, the 2^ ~^ term is only a moderate constant which does not scale with the 
size of the network in most of the interesting scenarios. 



J2 min(d^,|§;„n(Ar^(H™„))^|) • 



mm{\Sin^^{U,nn)^\,dn) 



(n,m)£G 




(11) 



Cartesian product. 
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J\f(H.m'n), where m' is an arbitrary element in Kgub- Update C = C + 1. Repeat 
this process until M{{iimn '■ G IKs„fe}) = {0}, 'i^sub ^ IK^ with \^sub\ = C or 

c = \KV 

- For every 'Ksub ^ with A/'({H^„ : m G Ks„4) 7^ {0}, set A/'({H^„, : m G 
U ({1, 2, K}\K^)}) = M{{Yi^n : m G 
For every Rx-m, m G {1, 2, /sT}, calculate A/'^({H^„ : n G K,„4) = (n„eK,„X'^(Hr 
(A/'^(H„„))-L, K^^fe C {1, 2, using a similar process. 

• Step3 Design Subspace constraints §^ and SJ!„: For every Tx n, n G {l,2,...,Er}, 
generate a series of potential subspace constraints G {d„,(i„ + 1, A^*} with 
|S^(d)| = d, based on the principe that a subspace which has higher null space "weight'^ 
is selected with higher priority. Choose the subspace constraint §^ from the potential 
subspace constraints: = S\^{d*), where: 

d* = arg max dn(d — dn) — 

d&{d„,dn + l,-,Nt-W{'H.nn)\} 

{1,2,...,K} 

J2 minK, |(Ar^(H„„))^|)min(|§^,(d) n Ar(H^„)^|,dO- (12) 

For every Rx m, m G {1, 2, i^}, use a similar process to generate E>l^^{d), d G 
{dm, dm + I, Nr} and set SJ„ = where 

d* = arg max dm(d — dm) — 

mm 

)l} 

{1,2,...,K} 

J2 min(d^, l^l^id) n (Ar^(H„„))^|) min(|§^ n (Ar(H^„))^|, rf„)-(13) 

• Step4 Low complexity feasibility checking: Denote vl^,Vm, n,m G {1,2,...,K} as 
the number of the freedoms at Tx n and Rx m, respectively. Set = — 

dn), = dm{\^m\ — dm)- Dcnotc Cmn, u ^ m E {1,2, K} as the number of 
constraints required to eliminate the interference from the Tx n to the Rx m. Set Cmn = 
min {dm, n (Ar^(H„„))^|) min (d^, |§*, H A/'(H„„)^|) and w = 0, Vn ^ m G 
{1,2,..., K}. Use freedom-constraint assignment to check if the system is proper (Please 
refer to Appendix O for details.) If the network is not proper, go to Step 5. Otherwise, 
let W = {dl, d;, d*^} = D, = = S'^, Mn, m G {1, 2, K}, and exit the 

algorithm. 

^The weight of A/'({Hmn : m G Ka^b}) is X^msK ^ From the left hand side of (lilt , this weight is the maximum 
number of lA constraints that one can mitigate by selecting a one dimensional subspace in A/'({Hm„ : m G Ks„(,}). 
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• StepS : Update D = {di, d2, ...da — l,..., dx} and go back to Step 3, where n is given 
by 

n = arg max ^(c^n + c„^ - 4„ - c'^J - {v^^ + < - t'*' - <') (14) 

ne{l,2,...,X} \ ^ / 
\m=l / 

where {f^', t)^'} and {c^„, c'„^}, m G {1, 2, K} are the number of freedoms and lA 
constraints under updated subspace constraints §*/ and SJi' given by (fT2l) and (fT3l) with 
(i^ = (i„ — 1, respectively. 
Remark 3.3 (stream assignment and subspace design criterion in Stage I Algorithm): As re- 
vealed in [6J, the lA feasibility condition is the major limitation of the DoF performance 
achieved by MIMO interference networks. Hence, in order to enhance the network DoF 
performance, subspace constraints {S^*,SJ^} and stream assignment D* are designed to 
alleviate the lA feasibility condition as much as possible. 

• Recall (l25l) . (ITTl) . we can see that (fT2l) . (fT3l) choose the dimension of subspace con- 
straints d* to maximize the difference between the number of freedoms in precoder (or 
decorrelator) design minus the number of lA constraints endued by the Tx (or Rx). 

• The right hand side of (fT4l) represents the number of lA constraint - freedom in transceiver 
design saved by removing one stream from Tx-Rx pair n. The higher this number, the 
more "constraint demanding" this stream is. Hence, it should be removed first so that 
the network can become lA feasible easier. ■ 

Lemma 1 (Property of the low complexity lA feasibility checking): A partially connected 
MIMO interference network with stream assignment D and potential signal subspace {S^, SJ^} 
policy is proper (i.e. satisfies (HD) if and only if it can pass the low complexity lA feasibility 
checking in Appendix |Cj Moreover, the worst case complexity of the proposed checking 
scheme is 0{K^). 

Proof: Please refer to Appendix O for the proof. ■ 
From Lemma [H and Definition 13. 1[ we have the following theorem: 



Theorem 3.1 (Property of S^} and W): Under Assumption 12. 1[ the potential signal 
subspace {S^*,§]^}, together with the stream assignment pattern W from Stage [T] form a 
proper system. 

2) Stage II: Precoder and Decorrelator Determination: In this section, we shall elab- 
orate the stage II processing, which determines the precoders {V„} and the decorrelators 
{Urn} under the subspace constraints determined in stage I. Specifically, since the subspace 
constraint can be represented by the constrained precoder and decorrelator as in (|9), ([TOj, 
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the precoder and decorrelator design is given by the following optimization objective which 
minimizes the total interference leakage power in the network: 



minu;„,v^ Ef=M.>o f trace ((U'„S;:H^.S^:v;)^(U'^S;;:H^.Si*V;)) . 



The following algorithm is guaranteed to converge to a local optimum [|7]|. 

Stage 2 (Precoder and Decorrelator Determination): Given {D*, S^*, S^}, determined by 



• Stepl Initialization : Denote and as the structure matrices for decorrelators and 
precoders: Um = U'^S^, V„ = S^V^. Set S^* and S'^ to be the aggregation of the 
basis vectors in and S[^, respectively. Randomly generate V'„. 

• Step2 Minimize interference leakage at the receiver side: At each Rx m, such that 



where u^(rf) is the d-th row of U'^, fd[A] is the eigenvector corresponding to the d-th 
smallest eigenvalue of A, d E {1,2, d^}. 
• Step3 Minimize interference leakage at the transmitter side: At each Tx n such 
that dl > 0, update V^: V^{d) = //, ^(U'„S;:H^„SJ:)^(U'„S;:H^„SJ:) , 

where v^((i) is the d-th column of V^, d G {1, 2, d* }. 

Repeat Step 2 and 3 until V'„ and U'^ converges. Set V; = S^V; and U,*„ = U'^S;,, 
Vn,m G {1,2,...,K}. 

Remark 3.4 (Backward Compatibility of the Proposed Scheme): When the system is fully 
connected (i.e. A/'(Hm„) = A/'^(Hm„) = {0}, Vn,m G {1, 2, iT}), it is easy to check 
in Algorithm 1, S^* = C^*^^ and = C^^^''; this means that in fully connected quasi- 
static MIMO interference networks, the proposed scheme reduces to the conventional lA 
schemes proposed in [6] and [|71- When Nt > Nr, the algorithm shall first utilize the null 
spaces on the Tx side and design S^* to null off part of the interference, which is similar 
to jSl. However, given a general partial connectivity topology {A/'(Hmn), A^^(Hmn)}, the 
algorithm generalizes the conventional interference alignment by dynamically combining the 
interference alignment and interference nulling approaches. ■ 



In this section, we shall derive the analytical DoF performance achieved by the proposed 
scheme in a symmetrical partially connected MIMO interference network. In fact, although 
the algorithm itself applies to general typologies, analyzing such cases can be prohibitively 



Stage I: 




IV. Performance Analysis 
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complicated. Since there are too many parameters in the general partial connectivity pa- 
rameters {A/'(H,mn), A/'^(Hmn)}5 iTi,n E {1,2, K} in Section Hl-Al we shall focus on 
a symmetrical K-pair partially connected MIMO interference network in which the partial 
connectivity is induced by both the path loss and the local scattering. 

Definition 4.1 (Symmetrical Partially Connected MIMO Interference Channels): Consider 
a K-pair partially connected MIMO interference network with the following configuration. 
Each Tx has Nt antennas and each Rx has A^^ antennas. Each Tx-Rx claims d^^^ = df data 
streams, Vn G {1, 2, K}. The partial connectivity states (elaborated below) are expressed 
in terms of three key parameters L, Ei and E2, which characterize the connection density, 
the rank of the direct links and the rank of the cross links, respectively. Please refer to 
Appendix |E] for the details. 

(span(eArj(g))) , g G Ei if: m = n 
AfCHmn) = { {spm{e]\r^{q)))^ , q e E2{n — m) if: < \n — m\ < L or \n — m\ > K — L (15) 
^Ntxi otherwise 



(span((Afihr„ng))^))^,gGEi if: 



m = n 



„ . \ < In — ml < L or In — ml 

span((A^hJ^7(g))^) , g G E2(n - m) if: i- 1 ^^^^ 

> K-L 



£ixNr otherwise 

where H^^, Ar and eN,{q) are defined in i^, & and respectively. h^7(s) G C^'^''^ 
is the s-th column of H^^. and E2(n - m) are subsets of {1,2, Nt}, with |Ei| = Ei, 
|E2(n — m)| = £'2, Vn, m. ■ 
Theorem 4.1 (Performance of the Partially Connected K -pair MIMO Systems): The pro- 
posed algorithm could achieve W = {df,df,...,df} if 

/ E^ + mm{E^,Nr) mm{E^, N^) \ 

^ "^"^ [miniK - 1, 2L) + 2 ' min(ir - 1, 2L)f + 1 j ' ^''^ 

Proof: Please refer to Appendix |F] for the proof. ■ 
Remark 4.1 (Interpretation of the Results): Note that the total DoF of the system is given 
by Kdf, using (fTTl) . the system can achieve a total DoF up to: 

El + min(£;i, A^,) min(£;i, Nr) 



K 



max 



(18) 



min(is: - 1, 2L) + 2' min(K - 1, 2L)|^ + 1 
where the first term and second term in the "max" operation are contributed by restricting the 
precoders and decorrelators in the subspaces and E>^* obtained in the stage I algorithm. 
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In the following, we shall elaborate various insights regarding how the partial connectivity 
affects the gain of the system. 

• The gain due to partial connection: In a K-pak fully connected quasi-static MIMO 
interference channel, the system sum DoF is upper bounded by _ partial 
connectivity improves this bound in two aspects: 1) Gain due to path loss: As path 
loss limits the maximum number of Rxs that each Tx may interfere, the total DoF of 
the system can grow on 0{K); 2) Gain due to spatial correlation: When the spatial 
correlation in the cross link is strong (i.e. small E2), a ^ factor gain can be further 
observed. 

• Connection density versus system performance: For large K, the DoF of the system 
(fTSi) scales with ~ ^{j)' which shows that the network density is always a first order 
constraint on the system DoF. 

• Rank of the cross links versus system performance: The system sum DoF (fT8l) is a 
(non-strictly) decreasing function of E2, which means system sum DoF grows when the 
rank of the cross links decrease. When E2 = 0, the achievable DoF in (fTSl) is reduced 
to: K mm(Ei, Nr) = A'rank(H„„), which means that all Tx-Rx pairs are using all the 
dimensions of the direct link for transmission. 

• Rank of the direct links versus system performance: The system sum DoF (fTSl) is a 
(non-strictly) increasing function of Ei, which means that the system sum DoF increases 
when the rank of the direct links increase. High rank direct links help to increase system 
performance from two aspects: 1) They increase the DoF upper bound that Tx-Rx pairs 
may achieve. 2) They increase the maximum number of free variables in the precoders 
and the decorrelators. 

• Backward compatibility with previous results: When the network is fully connected, 
i.e. El = E2 = Nt and L > [f], the inequality in ([iT]) is reduced to: df < which 
is consistent with the results in [|6l. ■ 

V. Simulation Results 

In this section, we shall illustrate the performance of the proposed scheme by simulation. To 
better illustrate how physical parameters such as the path loss and the scattering environment 
affect system performance, we consider the following simulation setup based on a randomized 
MIMO interference channel. 
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Definition 5.1 (Randomized Partially Connected MIMO Interference Channels): We have 
32 Tx-Rx pairs distributed uniformly in a 10A;m x 10A;m square as illustrated in Fig. |7] Each 
node has 12 antennas^ Each Tx-Rx pair is trying to deliver 2 data streams. Each Tx is 
transmitting with power P. Denote Dmn as the distance between the Tx n and Rx m. The 
partial connectivity is contributed the following factors: 

• Path loss effect: If Dmn > L, we assume the channel from the Tx n to the Rx m is 
not connected (H^n = 0). 

• Local scattering effect: If Dmn < L, then due to local scattering, the angular domain 
channel states = {hmnip,(l)}:P,<l ^ {1, 2, 12}, has the following property: 
hmn{p,C[) = if g satisfies where S is the radius of the local scattering otherwise 

hmn{p,q) --CAr{0,l). 

As a result, the partial connectivity parameters for the randomized model is given by 
»r/TT N I span(e^,(g)) , q G Qmn if- Dmn ^ L 

[tlmn) = \ (ly) 

C otherwise. 



^rHfT. ^ J (span((A^hr('Z)n) , q ^ Qmn if: Dmn < L 
j-^ixi2 otherwise. 



(20) 



where Umi^ and e^tiq) are defined in ©, © and (jS]), respectively, h^^(s) G 
is the s-th column of H^^, and is the set of all the column indices q G {1,2, 12} 
that satisfies Q. Note that Qmn is a random set with randomness induced by the random 
positions of the Tx n and Rx m. ■ 
Remark 5.1 (Physical Meaning of the Parameters in Definition \5.1\l : There are two param- 
eters in Definition 15. 1[ L and S. As as illustrated in Fig. |7l L is the maximum distance that 
a Tx can interfere (e.g. the big circle centered at Tx 1 in the figure) and hence reflects 
the connection density of the network. S is the radius of the local scattering, from (|6]), if 
the direction of a beam from the Tx does not overlap with the local scattering area of a 
Rx, it cannot be received by the Rx. (e.g. The local scattering area for Rx 1 is the small 
circles centered at Rx 1. Beam 7, 8 cannot be received by Rx 1 as their direction does not 
overlap with this circle.) Hence, S controls the rank (spatial correlation level) of the non-zero 
channels matrices. Larger S corresponds to higher rank channel matrices. ■ 

^We choose relatively large number of Tx-Rx pairs and antennas so that we can have a smooth system performance 
variation w.r.t. to partial connectivity parameters of the network such as L and S. 
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The proposed interference mitigation scheme is compared with 5 reference baselines below: 

• Conventional interference alignment (Baseline 1): The system directly adapts the 
precoder-decorrelator iteration proposed in 

• Maximum rank signal subspace (Baseline 2): Each node selects a maximum rank 
subspace constraint, i.e. set d* to be Nt - |A/'(H„„)| and Nr - \J\f^iiimm)\ in (fT2l) 
and (fT3]) . respectively in Stage I, then uses Stage II to determine the precoders and 
decorrelators. 

• Minimum rank signal subspace (Baseline 3): Each node selects a minimum rank 
subspace constraint, i.e. set d* to be dn and d^ in (fT2l) and (fT3l) . respectively in Stage 
I, in the stage I algorithm, then use the stage II algorithm to determine the precoders 
and decorrelators. 

• TDMA (Baseline 4) refers to the case where the Tx-Rx pairs use time division multiple 
access to avoid all interference. 

• Isotropic transmission (Baseline 5) refers to the case where the Tx and Rx sends and 
receives the data streams with random precoders and decorrelators without regard of the 
channel information. 

A. Performance w.r.t. SNR 

Fig. [8] illustrates the throughput per Tx-Rx pair versus SNR (lOlog^o(-P))- Here L = 5km 
and S = 3km. Conventional interference alignment (BL 1) saturates in the high SNR region 
as traditional lA is infeasible in this dense network. Both the proposed scheme and the 
Maximum/Minimum signal subspace methods (BL 2 and 3) can achieve throughputs that 
grow linearly with SNR since the on/off selection in the stage I algorithm guarantees that the 
system is feasible for lA. However, the proposed scheme achieves much higher DoF (^ 57) 
than BL 2 (^ 44) and 3 (~ 46), illustrating the importance of carefully designing the signal 
subspaces. Comparison of the proposed scheme and BL 1 shows that introducing subspace 
constraints can indeed enlarge the lA feasible region and enhance the system performance in 
both DoF and throughput sense. Moreover, note that for a 12 x 12, 2 stream per Tx-Rx pair 
and fully connected interference network, at most total network 22 DoF can be achieved, the 
performance of the proposed scheme (50 DoF) show that partial connectivity can indeed be 
exploit to significantly increase network total DoF. 
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B. Performance w.rt. Partial Connectivity Factors 

To better illustrate how different partial connectivity factors such as path loss and spatial 
correlation affect system performance, we illustrate the sum throughput versus L (the maxi- 
mum distance that a Tx can interfere a Rx) and S (the radius of the local scattering) under a 
fixed SNR (40dB) in Fig. [9] and Fig. [lOl respectively. By comparing the performance of the 
proposed scheme with different partial connectivity parameters, we have that the performance 
of the proposed scheme scales O {jjg), which illustrates a consistent observation as in Re- 
mark |4T|that weaker partial connectivity can indeed contribute to higher system performance. 
Moreover, comparison of the proposed algorithm with Baseline 2 and 3 illustrates how we 
should select signal subspaces and S[^* under different partial connectivity regions. For 
example, low rank subspace is more effective at high spatial correlation (small S) while high 
rank subspace is more effective at low spatial correlation (large S). Low rank subspace is 
also more effective compare to high rank subspace in dense networks (large L) and vice 
versa. By dynamically selecting signal subspace according to the partial connectivity state of 
the network, the proposed scheme obtains significant performance gain over a wide range of 
partial connectivity levels. 

VI. Conclusion 

In this paper, we have investigated how the partial connection can be utilized to benefit 
the system performance in MIMO interference networks. We considers a general partial 
connection model which embraces various practical situations such as path loss effects and 
spatial correlations. We proposed a novel two-stage interference mitigation scheme. The stage 
I algorithm determines the stream assignment and the subspace constraints for the precoders 
and decorrelators based on the partial connectivity state. The stage II algorithm determines the 
precoders and decorrelators based on the stream assignment and the subspace constraints as 
well as the local channel state information. The signal spaces is designed to mitigate "many" 
lA constraints at a cost of only a "few" free-variables in precoders and decorrelator so as 
to extend the feasibility region of the lA scheme. Analysis shows the proposed algorithm 
can significantly increase system DoF in symmetric partially connected MIMO interference 
networks. We also compare the performance of the proposed scheme with various baselines 
and show via simulations that the proposed algorithms could achieve significant gain in 
system performance of randomly connected interference networks. 
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Appendix A 

Physical Interpretation for Virtual Angle Model in MIMO Channel 

As has been observed in many previous works ffT4l . [fr5l the statistical property of the chan- 
nel states in a MIMO system is strongly affected by the physical propagation environment. 
For instance, in cellular MIMO systems, the Txs are positioned at high elevations above the 
scatterers while the Rxs are positioned at low altitude with rich scattering (Fig. [TK). Hence, 
only the scattering objects surrounding a Rx could effectively reflect signals from the Tx to 
the Rx as illustrated in Fig. [l^. Following an approach similar to [[T2l . [[T3l . we assume the 
double directional channel response from the n-th Tx to the m-th Rx (n, m E {1,2,..., K}): 
Hm„ = {h^mS^t.Or), 0t,9r G [-7r,7r)} has the following property: 

hlni0t,9r) = if \9t\>^; (21) 



2arcsin(;7^) when: S < dmn 
where a = { (22) 



, s 

27r else. 

where dmn is the distance between Tx n and to Rx m, S is the local ejfective scattering 
radius. Assume the Tx are equipped with uniform linear antenna array (ULA). Hence the 
virtual angular channel representation IfTTTl is given by: 

hmniP, Q) = J J hmniOt, dr)fNr ( W ) [ A iV* / ^^'"^^^'^^^^ 

where r is the antenna separation, A is the wavelength, 9mn is the angle between the transmit 
array normal direction and the direction from Tx n to Rx m, Lpmn is the angle between 
receive array normal direction and the direction from Tx n to Rx m, assume the antenna 
array is critically spaced, i.e. I = i. Mu) = eUO)e^iu) = ie-^---(^-^)^g=^. 

As illustrated in Fig. |2]\, fNi^) represents the radiation pattern of ULA [TS]. Note that 
the main-lobes dominate in the radiation pattern (e.g. when = 8, the power of the main- 
lobes occupy 91% of that in the whole radiation pattern). Hence, as illustrated in Fig. [2^, 
for simplicity, we use the main-lobe to approximate the radiation pattem, i.e. we use 



/7v(w), if: - [wj < ^ or > 



7V-1 



f'M={ --'-^ - ^ (24) 

0, otherwise. 

to replace fN{^) in (l23l) . Hence, combining (|2TI) . (I22l) and (1231) . we have 
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Appendix B 

The Number of Freedoms and IA Constraints in Equations Q and (fTOl) 
The freedoms in and U'„ of ^ and (flOl) are given by: 

dim{g{dn, 5"*)) = rf„(S'* - d/), and dim{g{dm, SU) = dmiS^ - d^), respectively. (25) 
where the Grassmannian Q{x, y) [fT6ll . IfTTll denotes the set of all x-dimensional subspaces in 

Then consider the number of independent constraints in (fTO^ . Consider the singular value 
decomposition of H„„ = U„„diag(si, Sa, sat) where U„„,V„„ are A^^ x A^^ and 
Nt X Nt unitary matrices, respectively, = mm{Nt, Nr), si ~ sn are the singular values 
of Hmn in descending order. Suppose rank(Hmn) = r, then we have: Hmn = [U^„, ^ 
diag(si, S2, Sr) 





respectively. Note that 



[V^„, V^^]-^, where and are A^^ x r and A'i x r matrices. 



U:„S;,H^.S*„V; = O ^ U:„S:;,ULdiag(.i, s^, Sr){VWSiV'^ = (26) 

andspan(VD = (span(VL))^ = (Ar(H„„))^, span^(UD = (span^(UL))^ = (Ar^(H^n))^, 
where span(X), span^(X) denote the linear space spanned by the columns of X and the 
rows of X^, respectively. Hence, the number of independent constraints in (flOl) is given by: 

rank(U:„S:;,UD ■ rank(VLS^V;) = min [dm, l^';^ n (Ar^(H^„))^|) min n Ar(H„„)^U„)(27) 

Appendix C 

Low Complexity Feasibility Checking Algorithm in Step 4 of Stage [H 

• Initialize the constraint assignment: Randomly generalize a constraint assignment pol- 
icy, i.e. {4„, d^J such that: 4^, c'^^ e NU{0}, + = c^„, m,ne {1, 2, 
(Note that in the subscripts of {c^^}' transmitter indexes come first). Calculate the 
variable - assigned constraint pressure, i.e. {P^,P^}, where 

me{l,2,...,i<'} nG{l,2,...,A'} 

• Update the constraint assignment: While there exist "overloaded nodes", i.e. -P* < 
or P^ < 0, m, n G {1,2, K}, do the following to update constraint assignment 

ic* c"" V 



January 20, 2013 



DRAFT 



21 



- A. Initialization: Select an "overloaded node" with negative pressure, without losing 
generality, assume this node is Tx-n, < 0. Set to be the root node of the 
"pressure transfer tree", which is variation of the tree data structure, with its nodes 
storing the pressures at the Txs and Rxs, its link strengths storing the maximum 
number of constraints that can be reallocated between the parent nodes and the 
child nodes. Please refer to Fig. |6]for an example. 

- B. Add Leaf nodes to the pressure transfer tree: 

For every leaf nodes (i.e. nodes without child nodes) {x E {t,r},n E {1,2,..., K}) 
with depths equal to the height of the tree (i.e. the nodes at the bottom in Fig. (6]): 

For every m G {1, 2, K}: If c^^ > 0, add P';^ as a child node of P';^ with link 

strength c^^, where x is the element in {t, r} other than x. 

- C. Transfer pressure from root to leaf nodes: For every leaf node just added 
to the tree in Step B with positive pressure, transfer pressure from root to these 
leafs by updating the constraint assignment policy {c^n,c^„}- For instance, as 
illustrated in Fig. |6^, F^^ "^""^ "^"^ P^^ is a root-to-leaf branch of the tree 
(red lines). Transfer pressure from P^^ to P^^ by updating: (c^^^J' = c^^^^ — e, 

(Cmim)' = ^mim + (Cmma)' ^ '^mma ~ (^nami)' = C^2mi + ^- HcnCC WC haVC 

(Pni)' = P-iii ~ ^ ^iid (-^4)' ^ -^4 + ^' where e is the minimum of the absolute 
value of the root pressure, leaf pressure, and all the strengths of the links, i.e. 
e = min (— , F^^ , c^^^^ , c^^nj , ^' denotes the value of A after update. Similarly, 
this operation can also be done for the green lines in Fig. |6B- 

- D. Remove the "depleted" links and "neutralized" roots: 

* If the strength of a link become after Step C: Separate the subtree rooted from 
the child node of this link from the original pressure transfer tree. 

* If the root of a pressure transfer tree (including the subtrees just separated from 
the original tree) is nonnegative, remove the root and hence the subtrees rooted 
from each child node of the root become new trees. Repeat this process until all 
roots are negative. For each newly generated pressure transfer tree, repeat Steps 
B~D (Please refer to Fig. |6]C! for an example). 

- E. Exit Conditions: Repeat Steps A~D until all trees become empty (hence the 
network is lA feasible) or no new leaf node can added for any of the non-empty 
trees in Step B (hence the network is lA infeasible). Exit the algorithm. 
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Appendix D 
Proof for Lemma [H 

We shall first prove the "if side. From Step 4 in Stage [T] and the Initialization step in 
Appendix O (fTTI) can be rewritten as: 

Yl ("^-rn + Cl^n) = < ^ + V^, (29) 

(n,m)GG (n,m)GG nSGn mSGm 

VG„, Gm C {1,2,..., K}, G = Gm ©G„. From the exit condition (Step E) of low complexity 
lA feasibility checking algorithm, we have > 0, > 0, Vn,n G {1,2,...,K}. Hence 
we have: 

neGn meGm (n,m)6G neG 

neG„ ms{i, meGm "6{1, 

2,...,K} 2,...,K} 

= Y.Pn+ P'rn>^ (30) 

neGn meGm 

VG„, Gm ^ {1,2,..., iT}, G = Gm © G„. This completes the "if side proof. 

Then we tum to the "only if side. We shall try to prove the converse-negative proposition 
of the original statement. If the network cannot pass the low complexity lA feasibility test, 
from the exit condition (Step E), there must exists a non-empty pressure transfer tree such 
that: 

• Root node has negative pressure. 

• All other nodes are non-positive. This is because positive nodes are either "neutralized" 
by the root in Step C if the strength of the links from the root to these nodes are sufficient 
or separated from the tree in Step D if one of the link strength is not sufficient. 

• No other nodes can be added to the tree, which implies = and c^,„, = for any 
Tx-n, Rx-m' in the tree and Rx-m, Tx-n' not in the tree. 

Hence, set G„, Gm in (l29l) to be the indexes of the Txs and Rxs that are in the remaining 
pressure transfer tree, we have: 

^ ^ (^n ^ ^ ^nm) ~^ ^ ] (^m ^ ^ ^mn) ^ ] (^n ^ ^ '-'nm) ~^ 

neGn meGm meGm neGn neGn inG{l,2,...,K} 

E o=E^"+E^™<o (31) 

meGm ne{l,2,...,X} neGn meGm 

Hence, the network does not satisfy (|29|) . This completes the "only if side proof. 



January 20, 2013 



DRAFT 



23 



Finally, let us consider the complexity of the checking algorithm. Since there are only 2K 
nodes, the algorithm can at most generate 2K trees. For each tree, since each cross link can 
be added into the tree once, there are at most K{K — 1) times of adding node operation. 
Hence the worst case complexity is 0{K^). 

Appendix E 

Detail Modeling of Symmetric Partially Connected MIMO Interference 

Network in Definition 14.11 

• Partial Connectivity due to Path Loss: If L < |n — m| < K — L assume H^n = 0. 

• Partial Connectivity due to Local Scattering (Direct Link): If n = m, assume due to 
local scattering, h^7(p) = if: j9 ^ Ei, otherwise h^7(p) ~ CN^''^\ where h^7(p) 
is the p-ih. column of the angular representation of the channel state HJ^^ (defined in 
©), p e {0, 1, Nt - 1}, El C {0, 1, ...,Nt-l} are the indices of the "good" angles. 
Denote |Ei| = Ei. 

• Partial Connectivity due to Local Scattering (Cross Link): If < |n — m| < L or 

\n — m\ > K — L and n ^ m, assume due to local scattering: h^^(p) = if: p ^ E^", 
otherwise h^^(p) ~ CA/"^""^*, where An = n - m e {-L,-L + 1,...L}, E^" are 
random subsets of {0, 1, A^* - 1} which satisfy |Ef "| = E2, VAn. < E2 < Nt. 

Appendix F 
Proof for Theorem 14. II 

Due to the symmetry property of the system, when W = {dj, dj, df}: Ek* = S^, 
^n* = ^m' ^'^i ^ {I5 2, K} and hence the system satisfies (dB if and only if: 

K K 

J2 Yl (1^- ^ (Ar^(H„„))^U/) min n AT {Hmn)^\, df) 

n=l mj^n 

K K 

<Y^df{s':-df) + Y,dAsz-df) 

n=l m=l 

K 

^ min [S''* - |§"* n A/'^(H^i)|, c//) min (5** - |§** n A/'(H^i)|, c//) 

m=2 

<df{S'* + S'*-2df) (32) 

where §** = = Vn G {1,2, ...,K}, S'* = |§**|, S'* = 

In general, due to the randomness in A/'^(Hmi) and J\f{iimi) it is hard to obtain optimal 
S'** and S'^*. To obtain a fundamental insight, we shall consider two extreme policies: S*** = 
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df,S^* = mm{Nr, El) (smallest subspace dimension on the Tx side and largest subspace 
dimension on the Rx side) and S^* = Ei,S^'* = mm{Nr, Ei) (largest subspace dimension 
on both the Tx and the Rx side). On the receiver side, when S"^* = min(A^r,^i), S^* = 
{Af" (Hu))-^ . On the transmitter side, from the §^((i) obtained in Step 3A, Algorithm!!] 
shall have the following form: S^(d) = span (^ejvj(^), e7Vt(^), ...ejy^^^J^, where pd is the 
p-th index in P = {0, 1, A^t — 1}, in which the elements are ordered w.r.t to the metric 
Em=2 H^Nti^^) e A/'(Hi„)) in descending order. When S^* = dj, ^ become: 

K 

J2^f- |S" n Ar^(H^i)| < min(Ar,,,^i) - df 

m=2 

K 

^ df{K - 1) - X] ^ l(ejv,(^) G A/'(Hi„)) < min(iV,, E^) - df 

d=l m=2 * 
, Nt K 

<= ^/(^ - 1) - ]^ E E ) ^ -^(Hin.)) < min(iV., E,) - dj 

* d=l m=2 * 

<= df (^m{K - 1, 2L)^^ < min(Ei, A^^) - df (33) 
When d* = El, ^ is simplified to: 

min(ir - 1, 2L) ■ min(^2, df) < Ei + min(Ei, A^^) - 2df 
<= mm{K -1,2L) ■ df < Ei + mm{Ei,Nr) -2df (34) 
From (|33T ) and (|34] ). (fTTT ) is obtained, which completes the proof. 
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1 (-1) 1 (-1) 

A. Original Ul-A radiation pattern B. ULA radiation pattern after approximation 

(with both main-lobes and side-lobes) (with main-lobes only) 



Fig. 2. Radiation pattern of ULA before and after approximation. 
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Fig. 3. Partial connectivity state (Tx side) of a 2 x 2, 5-pair interference channel. 



Rx-1 
(1) 

Tx-1 
(1) 


Rx-2 
(1) 

1 


Rx-3 
(1) 

1 


Rx.4 
(1) 

1 


Rx-5 Rx-1 
(1) (1) 
Tx-1 
■ ° 


Rx-2 
(1) 

1 


Rx-3 
(1) 

1 


Rx-4 
(1) 

■ 


Rx-5 Rx-1 
(1) (1) 
Tx-1 

■ (0) 


Rx-2 
(1) 

1 


Rx-3 
(1) 

1 


Rx-4 

(1) 




Rx-5 
0) 




Tx-2 
(1) ' 





1 


1 


Tx-2 





1 


1 


Tx-2 

1 10) 





1 


1 





Tx-3 , 
(1) 


1 





1 


^ y ^ 


1 





1 




1 





1 


1 


Tx-4 
(1) ^ 


1 


1 







1 


1 










1 





1 


Tx-5 
(1) 


1 


1 


1 


° ^,if ^ 


1 


1 


1 







1 









A B 



Fig. 4. The number of the freedoms in precoder and decorrelator design versus the number of the remaining constraints 
in the 2x2, 5-pair interference network before and after subspace design (Step 1^3 of the Algorithm in Section Ull-Bt . 
The numbers in red and blue denote the remaining freedoms in the corresponding decorrelators and precoders, respectively, 
and the numbers in purple denote the number of the remaining constraints to null off the interference. 
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A1: F-V pressure with random assignment 
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B1: F-V pressure witli adjusted assignment 
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Fig. 5. The freedom - variable pressure at each node (i.e. P,*, P^, n,m £ {1,2,3,4,5}, please refer to ( 128) for 
mathematical definition) before and after adjusting constraint assignments (Step 4 of the Algorithm in Section IIII-Bb . 
The overloaded nodes and the nodes with extra freedoms are marked out with yellow and green color, respectively. The 
reassignment processes are highlighted using red boxes. 




Fig. 6. Illustrative example of the "pressure transfer tree" and the corresponding operations in Appendix |C] A) A tree 
generated in Step A and B; B) Pressure transfer in Step C; C) Removal of depleted links and neutrialized roots in Step D. 
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Fig. 7. An illustration of a randomized partially connected MIMO interference network. In this illustration, K — 32, the 
Txs and Rxs uniformly distribute in a lOfcm x 10km square. L is the maximum distance that a Tx can interfere and S is 
the radius of the local scattering. 




Fig. 8. Throughput per Tx-Rx pair versus SNR for the proposed algorithm (and the 5 baselines) in a randomized partially 
connected MIMO interference channel. The parameters are given by L = 5km and S — 3km. 
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Throughput Per Tx-Rx Pair Versus L {Connection Density of the Networli) 

The Proposed Scheme 




3.5 4 4.5 5 5.5 6 6.5 

L (km) 



Fig. 9. Throughput per Tx-Rx pair versus L(The maximum distance that a Tx can interfere) for the proposed algorithm 
(and the 5 baselines) in a randomized partially connected MIMO interference channel. The transmit power is given by 
10logio{P) = iOdB and the local scattering radius S = 2.5km. 




Fig. 10. Throughput per Tx-Rx pair versus S(Radius of local scattering) for the proposed algorithm (and the 5 baselines) 
in a randomized partially connected MIMO interference channel. The transmit power is given by 10logiQ{P) — 40dB and 
L = 5km. 
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